Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  384 1600 (N27E65 L=205)  1518     99        5 ( 344/1558- 424/1648) -60.3 E71-E59 N23-N30
      2  747 1552 (N30E29 L=242)  1255    231        1 ( 705/1453- 815/1643)  90.0 E31-E26 N22-N37
      3 1667 1527 (N43W49 L=320)   839    114        0 (1621/1483-1692/1557) -15.9 W44-W53 N39-N45
      4  427 1394 (N15E47 L=223)   438    125        0 ( 413/1342- 446/1458) -60.3 E49-E44 N11-N18
      5 1047 1365 (N22E03 L=267)   523    223        0 (1001/1267-1072/1430) -55.0 E07-W00 N15-N25
      6  180 1090 (S07E65 L=205)   789    173        1 ( 154/1036- 203/1130)  56.3 E69-E61 S10-S04
      7 1785 1099 (N15W54 L=324)   353     57        0 (1759/1076-1808/1116) -11.3 W51-W56 N14-N16
      8  235 1050 (S09E58 L=213)  1160    127        4 ( 197/ 998- 258/1105)  81.9 E62-E55 S12-S05
      9  328 1005 (S10E47 L=223)   414     60        0 ( 299/ 984- 353/1029)   0.0 E50-E45 S11-S08
     10  366  965 (S12E44 L=227)   433    123        1 ( 322/ 934- 398/1023)  45.0 E47-E41 S14-S08
     11  255  943 (S15E55 L=215)   470     63        0 ( 231/ 911- 280/ 971) -21.8 E57-E51 S16-S13
     12  203  749 (S29E66 L=205)  1315    103        2 ( 157/ 714- 249/ 778)  12.3 E80-E57 S31-S26
     13 1018  640 (S23W05 L=275)   320     58        0 ( 997/ 612-1044/ 680) -36.9 W03-W07 S25-S20
     14 1793  615 (S14W69 L=339)  1545     85        0 (1762/ 570-1828/ 653) -63.4 W63-W75 S17-S11
     15  921  584 (S29E00 L=270)   217     73        0 ( 883/ 572- 958/ 598)   0.0 E02-W02 S29-S27
     16 1320  479 (S30W30 L=301)  1526    377        1 (1154/ 440-1494/ 527)  23.5 W19-W42 S35-S24
     17 1187  249 (S50W34 L=304)   481     69        1 (1145/ 235-1224/ 269)  14.0 W29-W37 S51-S48
     18 1378  230 (S48W57 L=328)   730     52        1 (1350/ 206-1403/ 260)  -8.1 W50-W65 S50-S45
  

    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