Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1074 1652 (N39E11 L=309)  1543    347        2 ( 925/1563-1204/1731) -13.5 E25-W00 N35-N42
      2  505 1581 (N22E49 L=271)  1356    158        2 ( 470/1502- 548/1653) -55.0 E55-E45 N17-N27
      3  247 1483 (N12E73 L=248)  3721    227        5 ( 190/1397- 301/1586) -69.0 E79-E66 N06-N19
      4 1577 1279 (N25W32 L=352)   531    128        0 (1533/1220-1614/1334) -45.0 W27-W34 N21-N27
      5 1153 1286 (N16W03 L=323)   320     76        1 (1122/1258-1174/1329)  90.0 W01-W04 N14-N19
      6  216 1081 (S12E59 L=262)  1020    184        3 ( 180/1016- 251/1135)  76.0 E63-E54 S16-S08
      7  280 1004 (S16E51 L=270)   918    248        1 ( 209/ 918- 364/1072)  38.7 E58-E43 S22-S10
      8  388 1002 (S14E41 L=279)   270     93        0 ( 369/ 963- 415/1043) -45.0 E42-E38 S15-S11
      9 1614  875 (N01W40 L=  0)   457     73        1 (1574/ 845-1650/ 901)  -8.1 W36-W43 N00-N02
     10  341  818 (S26E44 L=276)   394    112        0 ( 303/ 775- 386/ 865)  41.2 E48-E39 S29-S22
     11 1161  630 (S22W16 L=337)  1691    429        0 ( 985/ 544-1363/ 686)  35.5 W07-W28 S30-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