Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  912 1747 (N49W04 L=339)   276     62        0 ( 878/1732- 942/1759)   0.0 W01-W07 N47-N51
      2  347 1533 (N42E54 L=280)   475     42        0 ( 325/1510- 372/1553)   0.0 E60-E49 N39-N44
      3  582 1298 (N22E25 L=309)  1982    377        4 ( 442/1172- 691/1379) -26.6 E37-E18 N13-N29
      4 1524 1338 (N12W36 L= 11)   280     54        1 (1493/1318-1556/1360) -26.6 W35-W38 N10-N14
      5  143 1207 (N22E69 L=265)  1148    193        2 ( 101/1178- 186/1238)   7.6 E81-E61 N20-N25
      6 1696 1033 (S08W44 L= 18)   287     64        0 (1670/1011-1715/1065)  36.9 W42-W45 S10-S06
      7 1417 1010 (S06W24 L=358)   194     48        1 (1391/ 997-1449/1029) -26.6 W22-W25 S07-S05
      8 1875  957 (S15W62 L= 37)  1613    221        2 (1839/ 881-1916/1028)  38.7 W58-W70 S20-S10
      9 1345  969 (S08W18 L=353)   535    155        2 (1274/ 945-1414/1019) -32.0 W15-W23 S10-S04
     10 1218  865 (S12W10 L=344)   488    171        1 (1161/ 782-1262/ 953)  52.1 W07-W12 S18-S06
     11 1891  847 (S22W66 L= 41)   711    109        2 (1872/ 812-1913/ 896) -78.7 W63-W72 S24-S19
     12  212  668 (S11E67 L=267)  7043    765        8 ( 139/ 541- 357/ 791) -56.3 E73-E58 S20-S03
     13 1624  667 (S30W39 L= 14)  1753    281        4 (1546/ 586-1669/ 780)  71.6 W32-W46 S36-S22
     14  361  692 (S11E50 L=284)   354     72        1 ( 339/ 654- 374/ 727)  53.1 E51-E50 S14-S09
     15  558  559 (S22E39 L=295)  1546    328        3 ( 509/ 438- 599/ 697)  46.8 E38-E42 S30-S13
     16 1391  550 (S35W20 L=354)   416     89        1 (1371/ 506-1408/ 592) -80.5 W18-W22 S38-S31
  

    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