Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  962 1784 (N59W14 L= 48)  1349    291        3 ( 823/1739-1091/1822)   0.0 E03-W31 N56-N60
      2 1533 1660 (N36W57 L= 91)   392     55        0 (1507/1643-1563/1680)   0.0 W53-W62 N34-N36
      3 1370 1511 (N30W33 L= 67)   625    242        0 (1273/1435-1445/1559) -54.0 W27-W36 N23-N34
      4  348 1346 (N34E48 L=346)  3387    481        3 ( 263/1136- 466/1518)  61.6 E54-E39 N21-N45
      5 1169 1458 (N30W17 L= 51)   339    123        1 (1116/1429-1227/1494) -39.8 W14-W20 N27-N33
      6 1815 1269 (N06W62 L= 96)  3441    458       10 (1767/1079-1899/1429)  53.6 W51-W76 S03-N16
      7  175 1226 (N27E68 L=326)  2284    233        6 ( 126/1137- 230/1302) -21.0 E85-E58 N21-N31
      8 1331 1213 (N13W22 L= 56)   261    116        0 (1289/1169-1370/1244)  26.6 W18-W24 N10-N14
      9  747 1175 (N19E15 L= 19)   214    118        0 ( 709/1144- 796/1216) -50.2 E17-E12 N16-N22
     10 1550  822 (S14W30 L= 65)   877    191        0 (1466/ 766-1631/ 874)  14.0 W24-W36 S16-S12
     11  695  601 (S15E27 L=  7)   931    320        1 ( 596/ 514- 752/ 715)  38.2 E36-E22 S19-S08
     12  395  624 (S11E50 L=344)   812    109        2 ( 374/ 567- 423/ 677)  60.3 E52-E46 S14-S07
     13  518  614 (S12E40 L=354)   349     64        0 ( 485/ 601- 557/ 630) -21.8 E42-E36 S13-S11
     14  577  565 (S16E37 L=357)   286     69        1 ( 564/ 532- 591/ 604)  68.2 E38-E35 S18-S13
     15 1542  433 (S42W31 L= 66)   373     60        0 (1511/ 416-1572/ 454)   9.5 W28-W34 S43-S40
     16 1373  337 (S46W16 L= 50)   794    117        3 (1323/ 308-1429/ 369)   0.0 W11-W21 S49-S43
     17 1202  283 (S47E00 L= 34)   395     29        1 (1181/ 261-1226/ 308) -68.2 E01-W02 S49-S44
     18 1159  227 (S52E05 L= 29)  1118    108        2 (1111/ 202-1219/ 253) -18.4 E09-W00 S55-S49
     19 1179  163 (S62E06 L= 28)   697     67        2 (1146/ 148-1215/ 182) -35.5 E10-E01 S65-S58
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y