Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1534 1673 (N28W60 L=166)  1603    238        1 (1453/1584-1626/1720)   0.0 W46-W76 N26-N29
      2 1274 1456 (N25W28 L=134)   712    146        0 (1217/1402-1339/1493)   5.7 W22-W32 N22-N26
      3 1009 1287 (N22W06 L=112)   678    232        0 ( 974/1184-1070/1370)  48.4 W03-W11 N16-N25
      4  247 1272 (N38E56 L= 50)   721    168        0 ( 235/1207- 264/1357) -68.2 E58-E52 N34-N44
      5  713 1209 (N25E14 L= 92)   211     72        0 ( 686/1171- 733/1241) -68.2 E15-E13 N22-N27
      6 1161 1200 (N14W12 L=118)   239     47        2 (1130/1180-1192/1218)   0.0 W10-W14 N12-N15
      7 1898 1122 (S13W65 L=171)  2021    180        6 (1872/1052-1937/1204)  37.9 W60-W72 S16-S08
      8 1792  998 (S15W49 L=155)  1467    241        3 (1759/ 918-1831/1094)  60.9 W45-W53 S20-S09
      9 1739  876 (S21W42 L=148)   389     51        1 (1713/ 848-1771/ 897) -38.7 W40-W45 S22-S18
     10 1068  858 (S03E02 L=105)   206     84        0 (1036/ 842-1110/ 884) -63.4 E02-W00 S04-N00
     11 1676  820 (S22W36 L=142)   292     37        1 (1649/ 795-1700/ 841)   0.0 W33-W38 S23-S20
     12 1680  686 (S31W35 L=141)  1425    399        3 (1556/ 603-1783/ 793)  13.4 W23-W46 S33-S27
     13 1356  607 (S26W09 L=115)   811    371        1 (1249/ 501-1447/ 726)  27.9 W00-W17 S30-S20
     14  919  671 (S10E15 L= 91)   268    107        0 ( 872/ 656- 973/ 695) -45.0 E16-E11 S12-S07
     15 1629  555 (S39W30 L=136)   828    141        0 (1585/ 501-1671/ 617)  36.9 W25-W34 S41-S35
     16 1197  500 (S28E03 L=103)   287     60        0 (1181/ 474-1220/ 534) -90.0 E04-E02 S30-S25
     17  380  456 (S13E62 L= 44)  1473    121        2 ( 343/ 391- 411/ 501)  28.3 E70-E57 S17-S10
     18  592  424 (S19E45 L= 61)   374     78        0 ( 578/ 378- 603/ 468)  50.2 E47-E42 S21-S15
  

    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