Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1188 1843 (N58E20 L=219)   960     75        0 (1160/1804-1213/1882) -48.0 E26-E15 N52-N62
      2 1392 1601 (N43W09 L=248)   473     97        0 (1363/1559-1429/1636) -29.7 W06-W13 N40-N45
      3 1590 1603 (N51W30 L=269)   368     49        1 (1573/1576-1610/1626)  68.2 W27-W32 N48-N53
      4 1758 1478 (N47W51 L=290)   520     64        1 (1737/1442-1783/1502)  68.2 W47-W56 N43-N48
      5 1407 1331 (N24W16 L=254)   228     87        0 (1383/1301-1439/1365) -21.8 W13-W18 N23-N25
      6 1757 1188 (N26W46 L=285)  1945    357        1 (1655/1043-1854/1319) -27.6 W35-W59 N19-N31
      7  259  956 (S25E53 L=186)   966    196        2 ( 218/ 877- 315/1029)  73.3 E57-E47 S29-S19
      8 1233  674 (S18W22 L=260)   237     99        0 (1185/ 661-1281/ 690)  26.6 W18-W24 S19-S15
      9 1180  520 (S28W24 L=263)  2180    481        3 ( 966/ 450-1385/ 582)  39.6 W12-W35 S38-S19
     10 1458  485 (S23W45 L=284)   335     54        0 (1434/ 464-1485/ 507)  76.0 W43-W46 S24-S20
     11 1494  307 (S32W64 L=303)  1395    136        1 (1429/ 285-1566/ 331)   8.5 W56-W76 S34-S30
     12 1034  154 (S61W50 L=288)   888     73        0 (1015/ 119-1052/ 194) -20.9 W41-W62 S64-S56
  

    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