Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  844 1536 (N38E12 L=118)   756    229        1 ( 793/1475- 917/1630) -50.2 E16-E05 N33-N45
      2 1785 1474 (N27W65 L=194)  1051    148        2 (1755/1427-1810/1523)  80.5 W62-W67 N24-N30
      3 1134 1321 (N22W08 L=138)   680    128        1 (1073/1266-1192/1374) -45.0 W04-W11 N18-N26
      4  381 1338 (N25E46 L= 83)   683    139        2 ( 323/1317- 430/1360)  18.4 E51-E42 N23-N26
      5 1551 1247 (N15W37 L=166)   655    125        1 (1501/1205-1596/1294)  29.7 W32-W40 N13-N18
      6 1296  709 (S17W16 L=146)   393     87        1 (1265/ 670-1324/ 745) -78.7 W14-W18 S19-S14
      7  604  421 (S35E36 L= 94)  1542    503        1 ( 514/ 319- 715/ 568)  36.5 E50-E23 S43-S23
  

    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