Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  803 1637 (N40E25 L=296)   383    101        3 ( 768/1605- 845/1666) -14.0 E28-E20 N37-N41
      2 1190 1406 (N27W08 L=328)  2165    616        7 ( 927/1249-1394/1571) -25.2 E12-W22 N18-N36
      3 1541 1472 (N35W35 L=356)   769    247        3 (1440/1429-1630/1516) -11.3 W27-W42 N32-N37
      4  413 1470 (N23E50 L=271)   687    135        2 ( 378/1434- 450/1504) -24.0 E54-E45 N21-N25
      5 1847 1310 (N26W62 L= 23)   835    218        0 (1807/1282-1889/1328)   5.2 W57-W68 N23-N27
      6  897 1221 (N13E10 L=311)   351    182        1 ( 856/1158- 916/1309)  90.0 E11-E08 N08-N18
      7  745 1107 (N04E18 L=303)   349     98        2 ( 707/1077- 787/1152) -31.0 E20-E15 N02-N06
      8 1461 1081 (N09W27 L=348)   299    111        0 (1426/1037-1502/1127) -39.8 W24-W30 N06-N11
      9  410  781 (S19E40 L=280)   870    272        3 ( 325/ 730- 504/ 840)  26.6 E48-E32 S23-S14
     10 1752  724 (S11W56 L= 16)   798    263        2 (1678/ 691-1822/ 754)  -7.1 W48-W64 S13-S09
     11 1505  549 (S24W39 L=360)   371    106        0 (1448/ 525-1557/ 576)  29.7 W35-W42 S26-S22
  

    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