Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  671 1784 (N32E50 L=214)  2168    253        5 ( 587/1714- 761/1873)  90.0 E54-E46 N25-N41
      2 1489 1802 (N64W17 L=281)   614     54        0 (1469/1774-1507/1827)  50.2 W13-W23 N58-N70
      3  946 1706 (N32E25 L=239)   298     39        1 ( 920/1685- 968/1730)  56.3 E26-E22 N30-N34
      4  773 1689 (N27E35 L=229)   771    180        1 ( 715/1660- 841/1721)  35.5 E39-E32 N24-N30
      5 1520 1452 (N31W21 L=285)  1694    505        0 (1378/1281-1704/1612) -20.1 W07-W37 N26-N39
      6  471 1545 (N11E50 L=214)   691    118        4 ( 420/1514- 518/1585)  50.2 E53-E46 N08-N14
      7  395 1568 (N12E58 L=206)   435     98        1 ( 355/1548- 443/1583)  12.5 E62-E53 N10-N13
      8  583 1480 (N09E38 L=225)   766    150        0 ( 544/1414- 610/1542) -49.4 E43-E35 N05-N12
      9 1773 1342 (N33W45 L=309)   649    105        0 (1723/1312-1815/1369)   0.0 W39-W50 N31-N34
     10 1182 1177 (N05W05 L=269)   191     67        0 (1150/1154-1211/1205)  76.0 W03-W06 N03-N07
     11  487 1138 (S12E35 L=229)   565    128        1 ( 439/1089- 538/1179)  50.2 E38-E32 S15-S09
     12 1907 1116 (N24W63 L=327)   792    115        1 (1870/1062-1938/1167) -14.0 W57-W69 N21-N25
     13  142 1004 (S25E69 L=195)  1554    149        3 ( 119/ 931- 173/1066) -82.9 E74-E64 S29-S21
     14  340  948 (S26E45 L=219)   457    106        0 ( 302/ 909- 386/ 986)  35.5 E48-E40 S28-S23
     15 1184  900 (S10W12 L=276)   218     95        0 (1130/ 889-1226/ 911)  38.7 W09-W14 S12-S08
     16  526  792 (S32E27 L=237)  2306    514        4 ( 343/ 663- 733/ 865)  32.0 E45-E12 S44-S23
     17 1573  518 (S20W52 L=316)   823    222        0 (1472/ 465-1660/ 570)  54.0 W47-W56 S25-S14
     18 1282  537 (S28W30 L=294)   517    169        1 (1214/ 520-1365/ 554)  15.3 W25-W36 S30-S25
     19 1480  371 (S31W55 L=319)  2318    304        1 (1329/ 344-1609/ 405)  27.6 W43-W66 S36-S24
  

    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