Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1137 1472 (N25W15 L=168)  1423    350        2 ( 942/1392-1277/1550) -39.8 W03-W22 N18-N33
      2  359 1491 (N41E44 L=109)   719     99        1 ( 322/1454- 400/1521)  38.7 E47-E39 N38-N42
      3  470 1495 (N39E33 L=120)   559     94        1 ( 420/1470- 519/1518)   0.0 E37-E28 N37-N40
      4 1753 1448 (N14W59 L=213)   647     98        1 (1717/1409-1794/1484) -50.2 W56-W63 N10-N16
      5 1471 1358 (N13W33 L=186)   364    100        1 (1444/1319-1494/1406)  45.0 W30-W35 N10-N15
      6 1146 1352 (N18W13 L=166)   226     68        0 (1102/1326-1176/1369) -45.0 W10-W14 N16-N19
      7 1619 1168 (S00W39 L=193)   461     97        0 (1603/1111-1638/1220)  78.7 W38-W40 S03-N02
      8 1358  589 (S31W15 L=168)  1009    265        0 (1225/ 535-1486/ 635)   0.0 W05-W24 S32-S29
  

    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