Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1254 1650 (N39E06 L= 97)  2364    499       10 (1088/1526-1361/1763) -17.7 E21-W04 N31-N43
      2  734 1632 (N22E37 L= 66)   261     53        1 ( 712/1610- 767/1650)  71.6 E38-E35 N20-N23
      3  888 1347 (N09E16 L= 86)   582    228        2 ( 850/1235- 940/1444) -42.3 E21-E10 N03-N13
      4  974 1395 (N14E13 L= 89)   690    259        0 ( 911/1318-1048/1442)  76.0 E16-E08 N08-N16
      5  180 1279 (S11E68 L= 34)   904    128        2 ( 151/1242- 203/1307)  53.1 E72-E64 S13-S09
      6  232 1241 (S12E60 L= 42)   876    169        4 ( 205/1161- 262/1306) -83.7 E62-E57 S17-S08
      7 1693 1205 (N25W37 L=139)   631    154        1 (1661/1162-1725/1254) -38.7 W33-W39 N22-N26
      8 1695 1128 (N20W38 L=141)   447     99        1 (1652/1109-1725/1156)  45.0 W35-W40 N17-N22
      9  397 1105 (S17E41 L= 62)   326     87        1 ( 362/1071- 428/1138)  59.0 E43-E38 S19-S14
     10  167  816 (S39E59 L= 43)   547    123        1 ( 144/ 799- 191/ 843)  31.0 E61-E56 S40-S37
  

    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