Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1107 1864 (N55E09 L=342)   622     82        1 (1072/1838-1161/1878)  20.6 E12-E04 N52-N58
      2  450 1708 (N33E61 L=291)   515     47        1 ( 438/1676- 466/1739) -39.8 E64-E57 N30-N35
      3 1171 1512 (N26W03 L=355)   502    155        0 (1137/1445-1194/1577) -71.6 W00-W04 N21-N30
      4 1782 1249 (N19W50 L= 42)  1087    227        1 (1743/1146-1808/1355) -71.6 W47-W53 N13-N25
      5  387  951 (S17E42 L=310)   365     78        0 ( 370/ 911- 408/ 985)  63.4 E42-E39 S18-S14
      6  200  853 (S23E62 L=290)   488     87        2 ( 172/ 824- 237/ 888) -23.2 E65-E56 S25-S21
      7  440  839 (S23E37 L=314)   277     83        0 ( 402/ 816- 459/ 873)  38.7 E40-E35 S24-S20
      8  616  783 (S25E24 L=328)   456    146        0 ( 580/ 717- 640/ 851)  71.6 E26-E22 S29-S20
      9  469  587 (S39E38 L=314)   362    106        0 ( 443/ 548- 512/ 640)  37.9 E41-E32 S41-S34
     10 1112  209 (S60W31 L= 23)   622     56        2 (1086/ 178-1141/ 244) -31.0 W27-W37 S62-S56
  

    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