Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1766 (N36E30 L= 90)   698    125        1 ( 819/1732- 942/1792)   0.0 E34-E23 N34-N37
      2 1402 1573 (N35W12 L=132)   967    238        2 (1271/1537-1489/1597)  10.6 W03-W19 N31-N36
      3 1114 1353 (N13E02 L=118)   235    120        0 (1083/1306-1151/1400) -38.7 E04-W01 N11-N15
      4  503 1251 (S03E36 L= 83)   350    126        0 ( 480/1196- 518/1297) -80.5 E38-E34 S06-N00
      5  402 1162 (S10E42 L= 78)   271    106        0 ( 386/1113- 413/1216) -81.9 E42-E40 S13-S06
      6 1679 1166 (N16W38 L=157)   666    146        0 (1619/1110-1719/1213) -24.0 W32-W41 N13-N17
      7  199 1157 (S13E61 L= 59)   557     84        2 ( 163/1130- 243/1182)   0.0 E65-E55 S14-S11
      8  769 1065 (S10E15 L=104)   634    327        0 ( 648/1008- 878/1124)  19.7 E22-E08 S14-S05
      9  207  995 (S22E58 L= 61)   765    106        1 ( 163/ 956- 257/1037)  26.6 E63-E52 S24-S19
     10 1734  940 (N05W46 L=166)   826    164        1 (1713/ 855-1748/1033) -71.6 W43-W49 N00-N09
     11 1207  879 (S11W13 L=133)   296     73        0 (1181/ 843-1227/ 918) -53.1 W11-W14 S13-S09
     12  946  755 (S24W01 L=120)   285    110        0 ( 900/ 736- 997/ 790)  -9.5 E02-W04 S25-S22
     13 1742  723 (S07W54 L=173)   419     51        1 (1721/ 694-1762/ 754) -76.0 W52-W56 S08-S04
     14 1700  641 (S12W54 L=173)   733    159        1 (1666/ 588-1723/ 693) -24.0 W48-W58 S14-S09
     15 1834  639 (S08W69 L=188)   922     97        4 (1805/ 592-1859/ 685)  90.0 W66-W73 S10-S05
     16 1766  633 (S11W61 L=180)   817    100        2 (1728/ 589-1808/ 673)  56.3 W57-W63 S13-S07
     17 1639  454 (S24W60 L=180)  1769    206        2 (1545/ 402-1722/ 510)  33.7 W52-W65 S26-S18
     18 1490  300 (S37W61 L=181)   655     70        1 (1458/ 270-1512/ 339) -45.0 W57-W66 S39-S34
  

    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