Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  529 1616 (N23E51 L= 40)   371     71        1 ( 497/1600- 567/1632)   8.1 E54-E47 N21-N24
      2 1590 1469 (N39W32 L=122)   733    156        1 (1527/1423-1661/1503)  -8.7 W25-W38 N37-N40
      3 1216 1366 (N22W04 L= 95)   223     42        2 (1196/1342-1237/1393) -45.0 W02-W05 N20-N24
      4 1077 1332 (N17E04 L= 87)   710    239        4 ( 994/1269-1150/1377)   0.0 E08-W01 N14-N19
      5 1320 1306 (N21W12 L=103)   332     54        0 (1291/1283-1342/1337)  76.0 W10-W13 N19-N23
      6 1325 1233 (N17W14 L=104)   548    153        2 (1277/1188-1369/1259)  -9.5 W10-W16 N13-N18
      7  974  600 (S28W07 L= 97)  1465    438        2 ( 785/ 533-1149/ 698)  43.5 E04-W15 S36-S18
      8  830  358 (S50W04 L= 95)   443     69        0 ( 800/ 337- 878/ 381)   8.1 W01-W08 S51-S48
      9 1105  292 (S47W31 L=121)  1739    233        2 (1001/ 267-1214/ 337)  23.6 W21-W39 S50-S43
  

    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