Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  915 1776 (N60W13 L=217)   770     89        1 ( 881/1743- 955/1799) -33.7 W08-W18 N57-N62
      2  282 1325 (N36E56 L=148)   566     86        1 ( 243/1301- 314/1347)   0.0 E62-E52 N34-N37
      3  852 1203 (N20E08 L=196)   718    273        2 ( 747/1174- 965/1243) -12.1 E14-W00 N17-N23
      4 1742 1030 (S07W48 L=252)   373     73        1 (1714/1003-1762/1069) -59.0 W46-W50 S09-S04
      5  155 1028 (N18E70 L=134)   989    103        1 ( 116/ 999- 204/1054)  -6.7 E78-E61 N15-N19
      6  988  789 (S07E06 L=198)   278     49        2 ( 960/ 768-1014/ 808) -56.3 E07-E04 S08-S05
      7 1318  746 (S15W13 L=217)   302    111        0 (1282/ 699-1356/ 788)  33.7 W10-W16 S17-S13
      8  272  651 (S06E63 L=141)  2098    407        3 ( 214/ 555- 347/ 727) -48.0 E68-E57 S12-S01
      9  482  646 (S08E43 L=161)   361     64        1 ( 461/ 619- 507/ 681)  36.9 E44-E40 S09-S06
     10 1787  618 (S37W56 L=260)  1295    132        1 (1727/ 589-1848/ 654)  -3.0 W47-W66 S38-S34
     11  920  565 (S20E15 L=189)   545    206        1 ( 864/ 494- 988/ 641)  36.0 E20-E09 S23-S15
     12  689  612 (S13E29 L=175)   244     68        0 ( 657/ 587- 715/ 636) -45.0 E30-E26 S14-S11
     13 1218  504 (S30W03 L=207)  1236    274        0 (1140/ 407-1324/ 606)  33.7 E04-W11 S35-S24
     14 1560  545 (S35W30 L=234)   597    106        1 (1507/ 520-1603/ 582)   0.0 W25-W33 S36-S32
     15  674  533 (S18E32 L=172)   707    193        0 ( 605/ 505- 742/ 580) -12.5 E35-E26 S20-S13
     16  531  418 (S24E49 L=155)   335     49        0 ( 516/ 387- 550/ 446)  38.7 E51-E45 S26-S22
     17  689  326 (S33E41 L=163)  3044    323        5 ( 635/ 178- 780/ 444)  25.5 E72-E28 S47-S24
  

    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