Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1594 1682 (N25W68 L= 63)   486     60        0 (1554/1674-1628/1694) -16.7 W63-W73 N22-N26
      2 1529 1614 (N24W55 L= 50)   548     86        1 (1492/1583-1564/1641)   0.0 W49-W59 N22-N25
      3 1396 1564 (N26W42 L= 37)   503    132        1 (1369/1502-1436/1615)  66.8 W39-W44 N21-N28
      4 1358 1419 (N19W32 L= 27)   243     78        1 (1333/1381-1380/1443)  21.8 W28-W33 N17-N20
      5 1141 1402 (N24W18 L= 13)   263     81        0 (1100/1384-1183/1422) -53.1 W15-W20 N22-N26
      6  993 1303 (N23W06 L=  1)   673    244        1 ( 881/1251-1079/1333) -38.7 E01-W09 N17-N25
      7  278 1205 (N35E51 L=304)   780    142        0 ( 225/1165- 344/1249) -21.8 E58-E43 N32-N38
      8  721 1171 (N23E14 L=341)   555    212        0 ( 692/1100- 772/1246) -84.3 E16-E12 N17-N27
      9 1912 1079 (S18W67 L= 62)  1330    129        1 (1853/1040-1957/1109) -21.8 W58-W78 S22-S14
     10 1641 1030 (S10W37 L= 32)   562    134        0 (1601/ 972-1693/1086)  26.6 W33-W42 S13-S08
     11 1850  993 (S20W56 L= 51)   741    127        3 (1824/ 939-1888/1040)  18.4 W52-W62 S23-S18
     12 1256  777 (S13W07 L=  2)   877    300        1 (1197/ 659-1315/ 913)  45.0 W00-W13 S18-S06
     13  500  836 (N09E36 L=319)   692    177        1 ( 454/ 775- 559/ 910) -71.6 E38-E33 N04-N13
     14 1573  814 (S21W27 L= 23)   285     64        1 (1551/ 784-1591/ 842) -68.2 W26-W28 S23-S18
     15 1272  522 (S30W01 L=356)   440     67        1 (1249/ 485-1299/ 558)  63.4 E01-W02 S31-S26
     16  821  425 (S23E30 L=325)   524    176        2 ( 778/ 356- 859/ 486)  35.0 E36-E26 S26-S19
     17  806  230 (S37E45 L=310)   536     81        2 ( 780/ 192- 836/ 271)  53.1 E48-E42 S41-S33
     18  983  224 (S43E33 L=322)   368     60        1 ( 963/ 200-1006/ 252)  23.2 E36-E29 S45-S41
  

    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