Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  756 1661 (N51W01 L=264)   441     71        1 ( 723/1629- 782/1697)   9.5 E01-W05 N48-N53
      2 1431 1656 (N29W49 L=312)   596     55        1 (1402/1625-1459/1681)   7.1 W45-W53 N26-N30
      3  277 1183 (N31E50 L=212)   679    114        0 ( 225/1142- 322/1218)  10.3 E56-E44 N29-N33
      4  652 1114 (N19E20 L=243)   632    208        0 ( 581/1076- 737/1161) -37.9 E23-E14 N15-N23
      5  403 1080 (N22E39 L=224)   340     92        2 ( 364/1059- 452/1105) -29.7 E41-E34 N20-N24
      6  304  952 (N17E49 L=213)   283     53        0 ( 286/ 928- 322/ 978) -71.6 E51-E47 N15-N18
      7 1343  729 (S20W11 L=273)   231    109        0 (1321/ 683-1376/ 766)  31.0 W08-W13 S22-S18
      8 1657  689 (S33W34 L=296)   291     31        0 (1632/ 674-1682/ 707) -11.3 W31-W36 S33-S31
      9  856  371 (S30E31 L=232)   303     57        1 ( 826/ 354- 886/ 400) -11.3 E33-E28 S31-S27
     10 1259  192 (S60E16 L=247)   706     68        0 (1224/ 174-1297/ 217) -71.6 E19-E12 S64-S55
     11 1175  164 (S60E27 L=236)   981     68        1 (1137/ 138-1216/ 185) -50.2 E33-E22 S63-S55
     12 1077  145 (S57E37 L=225)   570     31        1 (1055/ 129-1100/ 164)   7.1 E41-E32 S59-S54
  

    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