Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1123 1732 (N39E16 L= 62)   554     78        2 (1077/1708-1160/1751)   8.1 E19-E12 N37-N40
      2 1302 1663 (N39E01 L= 77)  1173    212        5 (1230/1611-1366/1714) -10.3 E06-W05 N35-N43
      3 1533 1622 (N45W19 L= 96)   371     53        1 (1509/1602-1561/1647) -11.3 W16-W21 N43-N46
      4  923 1524 (N19E19 L= 58)   673    190        0 ( 871/1454- 987/1604) -28.6 E24-E13 N16-N22
      5  251 1476 (N03E69 L=  8)  2519    375       10 ( 189/1426- 293/1554)  69.4 E77-E61 N00-N08
      6 1371 1457 (N27W09 L= 87)   581    146        2 (1336/1406-1399/1529) -45.0 W05-W12 N23-N31
      7 1411 1287 (N18W16 L= 93)   238     84        0 (1392/1239-1426/1327) -59.0 W13-W17 N14-N19
      8 1791 1202 (N25W45 L=123)  1439    256        0 (1708/1117-1861/1291) -16.4 W36-W53 N22-N27
      9 1570 1020 (N07W32 L=109)   344    149        0 (1519/ 974-1624/1051)   0.0 W28-W36 N04-N08
     10  263  924 (S30E50 L= 28)   632    148        2 ( 216/ 894- 312/ 952)   6.3 E54-E45 S31-S27
     11  747  444 (S50E01 L= 76)   640    137        1 ( 691/ 415- 824/ 485)   9.5 E06-W06 S51-S47
     12  859  213 (S64W25 L=103)   479     51        1 ( 835/ 196- 886/ 238)   8.1 W21-W29 S66-S62
  

    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