Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  462 1650 (N35E58 L=328)  1180    149        4 ( 394/1617- 519/1678)  -2.9 E68-E48 N33-N37
      2  559 1537 (N28E41 L=345)   671    128        2 ( 539/1462- 582/1595) -48.4 E45-E36 N23-N32
      3  372 1456 (N21E54 L=332)   745    191        0 ( 307/1397- 421/1544)  70.0 E58-E51 N16-N27
      4  317 1500 (N23E63 L=322)   521     65        1 ( 292/1476- 337/1526)  63.4 E66-E59 N21-N25
      5  906 1364 (N21E11 L= 15)   533    224        4 ( 864/1306- 982/1441) -33.7 E14-E05 N17-N25
      6  604 1068 (S00E27 L=359)   214     98        0 ( 560/1050- 646/1089)   0.0 E29-E23 S01-N00
      7 1082 1051 (N03W03 L= 29)   470    163        2 (1010/1021-1154/1081)  15.9 W00-W07 N01-N05
      8 1462  736 (S12W31 L= 57)  1645    493        2 (1302/ 624-1599/ 923)  53.1 W21-W38 S20-N00
      9  240  676 (S29E61 L=324)   682     90        0 ( 207/ 639- 278/ 708)  19.7 E68-E54 S31-S26
  

    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