Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1055 1611 (N44W10 L=130)   253     16        0 (1032/1596-1081/1626) -14.0 W08-W12 N42-N44
      2 1495 1584 (N35W46 L=166)   726    109        1 (1471/1542-1530/1621)  24.0 W41-W51 N32-N37
      3 1207 1386 (N26W16 L=137)   460    176        0 (1149/1338-1278/1447) -48.8 W13-W20 N22-N30
      4 1630 1304 (N15W46 L=166)  1270    280        5 (1569/1214-1686/1399) -90.0 W41-W50 N09-N21
      5  818 1103 (N12E12 L=108)   316    106        1 ( 781/1083- 870/1118)  -9.5 E14-E08 N10-N13
      6 1750  856 (S14W50 L=170)   488     38        0 (1729/ 830-1769/ 889)  76.0 W47-W51 S15-S11
      7  541  846 (S02E33 L= 87)   486     76        1 ( 514/ 809- 577/ 877)  36.9 E34-E29 S03-N00
      8  297  600 (S17E62 L= 59)   873     89        2 ( 264/ 563- 327/ 644) -71.6 E65-E57 S19-S13
      9  756  380 (S35E29 L= 92)   894    196        3 ( 705/ 316- 821/ 468) -76.0 E30-E25 S40-S27
     10 1189  449 (S34W05 L=126)   434    127        0 (1129/ 435-1251/ 464) -18.4 W01-W10 S35-S32
     11 1136  206 (S57E02 L=118)   760    114        1 (1084/ 188-1202/ 222)   0.0 E07-W05 S58-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