Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  665 1774 (N33E48 L=359)   711     76        2 ( 631/1745- 709/1800) -11.3 E52-E42 N31-N35
      2 1243 1756 (N46E04 L= 43)   562     88        1 (1212/1727-1279/1788)  74.1 E05-E01 N42-N49
      3  453 1667 (N22E59 L=349)  2644    202        4 ( 352/1618- 541/1736)  18.4 E70-E52 N18-N27
      4  881 1647 (N28E24 L= 23)   617    221        2 ( 782/1579- 959/1700) -14.0 E33-E17 N24-N30
      5 1511 1465 (N31W22 L= 69)  1961    492        2 (1322/1309-1648/1594) -20.3 W06-W34 N24-N34
      6 1126 1398 (N17E02 L= 45)   319    129        0 (1091/1349-1168/1442) -33.7 E05-W01 N14-N18
      7  667 1243 (S01E25 L= 22)   314     82        0 ( 638/1216- 687/1288)  90.0 E26-E24 S03-N01
      8 1825 1118 (N18W52 L=100)  1610    243        6 (1778/1059-1898/1219) -18.4 W46-W62 N14-N23
      9 1834  971 (N10W56 L=103)  1180    156        1 (1809/ 889-1853/1054) -78.7 W53-W58 N04-N14
     10  503  930 (S22E31 L= 16)   347     86        0 ( 492/ 880- 514/ 976) -71.6 E32-E30 S25-S19
     11  165  840 (S32E67 L=340)  1280    198        4 ( 131/ 743- 214/ 917)  26.6 E81-E58 S38-S27
     12 1709  647 (S11W56 L=103)   378     46        1 (1691/ 624-1729/ 675) -71.6 W53-W58 S12-S09
     13  524  628 (S41E28 L= 20)   283     74        0 ( 490/ 608- 564/ 642)  23.2 E30-E23 S43-S39
     14 1430  517 (S26W40 L= 88)  1050    282        2 (1332/ 474-1555/ 575)   2.6 W30-W52 S28-S23
     15 1459  335 (S36W56 L=103)  1384    139        2 (1409/ 292-1518/ 407)  -3.0 W45-W64 S38-S32
     16 1379  278 (S42W56 L=103)   554     56        0 (1359/ 248-1401/ 313) -16.7 W51-W61 S42-S39
  

    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