Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  656 1231 (N26E17 L=146)   859    204        2 ( 630/1136- 681/1335)  80.5 E19-E14 N20-N32
      2  807 1230 (N22E07 L=156)   266     86        0 ( 768/1203- 852/1246)   0.0 E10-E04 N20-N22
      3  159 1040 (N25E65 L= 98)   755     95        5 ( 133/ 998- 177/1091)  63.4 E69-E62 N22-N28
      4  197 1060 (N26E60 L=104)   382     86        1 ( 183/1013- 226/1091) -90.0 E61-E55 N23-N27
      5 1574  853 (S20W29 L=192)   348    122        1 (1526/ 805-1624/ 889)  14.0 W24-W32 S21-S18
      6 1706  646 (S37W38 L=201)   401     34        1 (1683/ 623-1726/ 672)  45.0 W35-W40 S39-S35
      7 1465  605 (S32W16 L=180)   424    163        1 (1395/ 556-1528/ 640)   5.2 W10-W21 S32-S30
      8 1235  402 (S38E05 L=158)   486    128        0 (1185/ 372-1301/ 446)   6.3 E09-W00 S39-S37
      9  540  271 (S29E65 L= 99)  1060     74        3 ( 506/ 229- 563/ 309)  20.6 E74-E58 S32-S25
  

    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