Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1606 1624 (N28W63 L= 75)   420     50        0 (1579/1603-1631/1645) -63.4 W60-W64 N25-N29
      2 1666 1592 (N24W66 L= 78)   571     54        1 (1640/1569-1686/1614)  15.9 W63-W71 N22-N25
      3  618 1440 (N40E21 L=351)  1044    184        1 ( 550/1398- 713/1499) -35.5 E27-E13 N35-N45
      4 1282 1331 (N19W22 L= 35)   222     92        0 (1267/1285-1301/1370)  78.7 W21-W24 N16-N21
      5 1037 1331 (N25W07 L= 19)   199    123        0 (1000/1298-1090/1346)   0.0 W04-W10 N23-N25
      6  170 1286 (N36E73 L=299)  1535    131        3 ( 141/1224- 209/1341) -29.7 E86-E64 N31-N39
      7  293 1132 (N25E51 L=321)   533    124        2 ( 256/1084- 316/1186) -45.0 E55-E48 N22-N28
      8  225  991 (N17E59 L=313)   615    166        0 ( 181/ 930- 258/1025) -45.0 E64-E55 N13-N19
      9  314  934 (N13E50 L=322)   321     99        0 ( 293/ 886- 333/ 979) -71.6 E51-E48 N09-N15
     10  404  920 (N11E42 L=330)   347     73        3 ( 370/ 889- 447/ 939)   0.0 E45-E38 N09-N12
     11  166  884 (N11E69 L=303)   538     72        2 ( 142/ 852- 188/ 919)  45.0 E72-E64 N08-N13
     12  218  804 (N06E63 L=309)   508     58        0 ( 195/ 780- 243/ 827) -31.0 E66-E60 N04-N07
     13  425  540 (S13E52 L=320)  2312    503        5 ( 318/ 415- 528/ 666) -72.5 E58-E46 S22-S03
     14  610  628 (S09E34 L=338)   468     95        1 ( 578/ 591- 642/ 663)  45.0 E37-E31 S11-S07
     15  571  514 (S16E41 L=331)   810    236        2 ( 500/ 440- 663/ 579) -55.0 E44-E36 S22-S11
     16 1590  500 (S40W31 L= 43)   571     76        1 (1561/ 475-1617/ 529)   0.0 W28-W34 S42-S37
     17 1280  415 (S37W04 L= 16)   279     52        0 (1261/ 390-1304/ 449)  36.9 W02-W06 S39-S35
     18 1115  387 (S35E09 L=  3)   249     56        1 (1089/ 370-1141/ 409)   0.0 E10-E06 S36-S33
     19  769  314 (S34E36 L=336)   641    181        1 ( 708/ 259- 808/ 381)  26.6 E45-E29 S37-S29
  

    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