Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  695 1502 (N31E27 L=355)   529    180        0 ( 623/1438- 756/1530) -21.0 E33-E20 N27-N33
      2 1446 1469 (N33W30 L= 51)   331     60        1 (1423/1446-1476/1496) -45.0 W28-W31 N31-N34
      3 1063 1413 (N27W01 L= 22)   280    117        0 (1011/1371-1104/1445) -38.7 E02-W03 N24-N29
      4  899 1222 (N14E09 L= 13)   229     66        0 ( 868/1206- 938/1245) -14.0 E10-E06 N13-N15
      5  338 1161 (N07E48 L=333)   317     47        0 ( 315/1137- 361/1189)  63.4 E50-E46 N05-N09
      6 1625 1046 (N06W40 L= 62)   285    130        1 (1578/1006-1681/1067)  20.6 W36-W44 N03-N07
      7  389  930 (S07E42 L=339)   759    217        2 ( 362/ 838- 425/1020)  80.5 E44-E39 S12-N00
      8  273  815 (S15E54 L=327)  2792    434        6 ( 192/ 683- 316/ 985)  49.4 E69-E48 S24-S03
      9  822  816 (S11E12 L= 10)   239    115        0 ( 777/ 782- 868/ 851) -26.6 E14-E08 S13-S09
     10  651  695 (S20E23 L=358)  1107    254        4 ( 585/ 595- 713/ 799) -69.0 E28-E19 S26-S13
     11  718  736 (S17E19 L=  3)   371     58        1 ( 694/ 702- 748/ 765)  78.7 E20-E16 S19-S14
     12 1843  731 (S14W67 L= 89)   676     77        1 (1800/ 714-1880/ 746)   9.5 W61-W73 S14-S12
     13  318  577 (S31E57 L=325)  2204    354        4 ( 240/ 496- 451/ 680) -36.9 E68-E44 S35-S23
     14  898  538 (S29E07 L= 15)  3255    640        6 ( 639/ 433-1145/ 671)  25.9 E28-W09 S38-S19
     15  552  615 (S26E32 L=350)   497    120        1 ( 509/ 578- 612/ 643)  29.1 E36-E26 S28-S23
     16 1285  489 (S31W22 L= 43)   634    199        1 (1199/ 446-1370/ 531)  32.5 W15-W27 S34-S27
  

    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