Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1272 1648 (N39E04 L= 96)  2362    469        8 (1114/1503-1378/1756) -19.8 E19-W06 N31-N44
      2  483 1499 (N08E48 L= 53)   360     93        2 ( 442/1473- 513/1527)  53.1 E50-E45 N05-N09
      3  956 1356 (N11E13 L= 88)  1208    451        3 ( 880/1225-1074/1433) -42.3 E19-E06 N03-N17
      4  245 1247 (S12E59 L= 42)   962    223        3 ( 215/1156- 276/1336) -79.7 E62-E56 S17-S06
      5  193 1276 (S11E66 L= 34)   680    153        1 ( 165/1248- 222/1299)  45.0 E70-E62 S12-S09
      6 1720 1158 (N23W40 L=141)  1192    173        4 (1675/1097-1750/1235) -40.6 W36-W43 N17-N27
      7  259  980 (S27E50 L= 50)   422    101        1 ( 232/ 953- 294/1017)  38.7 E53-E47 S29-S24
      8  173  806 (S40E58 L= 42)   535    124        0 ( 145/ 787- 200/ 822)  26.6 E61-E55 S41-S38
  

    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