Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1461 (N15E18 L=238)  1206    294        1 ( 859/1337- 951/1600) -51.3 E25-E12 N07-N22
      2  355 1443 (N03E55 L=200)  1562    342        3 ( 233/1405- 460/1477)  13.4 E67-E46 N00-N05
      3 1513 1317 (N22W22 L=277)   360     90        1 (1494/1274-1524/1362) -63.4 W20-W23 N19-N25
      4 1047 1253 (N06E04 L=251)   608    155        1 ( 971/1201-1109/1312) -18.4 E09-W00 N04-N07
      5 1520 1022 (N05W29 L=284)  1393    294        2 (1473/ 864-1576/1158) -40.6 W22-W36 S02-N12
      6  229 1076 (S21E55 L=200)  1315    251        5 ( 161/1032- 269/1144) -22.6 E63-E50 S23-S16
      7  669  858 (S25E17 L=238)   810    276        1 ( 541/ 826- 787/ 908)  33.7 E25-E10 S29-S19
      8  273  782 (S38E49 L=207)  1489    419        1 ( 176/ 729- 394/ 817)  22.2 E59-E37 S42-S33
      9 1650  404 (S22W64 L=320)  1838    219        2 (1582/ 348-1698/ 475)  66.0 W59-W68 S26-S17
  

    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