Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1644 (N45W07 L=127)  1188    222        1 (1344/1563-1396/1729) -70.3 W04-W11 N39-N53
      2  464 1644 (N21E57 L= 62)  1123    196        2 ( 413/1594- 497/1714) -77.5 E60-E54 N17-N26
      3  892 1670 (N33E27 L= 92)   326     35        0 ( 872/1648- 910/1694)  90.0 E28-E25 N31-N35
      4  552 1554 (N17E45 L= 75)   359    139        1 ( 531/1500- 581/1619) -63.4 E46-E42 N14-N22
      5 1644 1582 (N51W36 L=156)   413     85        0 (1621/1560-1672/1600)  26.6 W32-W40 N48-N52
      6  955 1559 (N27E18 L=101)   265     57        0 ( 926/1545- 987/1577)  14.0 E20-E16 N25-N28
      7  639 1506 (N16E36 L= 83)   344    122        0 ( 612/1475- 675/1550) -26.6 E39-E33 N14-N18
      8  321 1495 (N09E63 L= 57)   492     93        2 ( 289/1463- 349/1531)  59.0 E65-E61 N06-N11
      9  717 1516 (N19E32 L= 88)   252     45        0 ( 692/1499- 746/1531)  56.3 E33-E29 N17-N20
     10  959 1376 (N16E12 L=107)   478    143        1 ( 944/1307- 977/1444) -58.0 E14-E09 N12-N20
     11  220 1206 (S10E61 L= 58)   964    158        4 ( 183/1153- 252/1255) -80.5 E65-E57 S12-S06
     12 1283 1221 (N14W10 L=130)   360     69        2 (1248/1193-1315/1252)  36.9 W08-W12 N12-N16
     13  289 1212 (S08E54 L= 66)   415     77        1 ( 265/1179- 319/1233) -26.6 E56-E51 S09-S06
     14 1114 1058 (N01W04 L=124)   242     48        2 (1086/1038-1139/1087)  63.4 W02-W05 N00-N02
     15 1733 1002 (N13W46 L=166)   497    107        1 (1705/ 956-1761/1046)  74.1 W44-W48 N09-N16
     16  627  766 (S27E19 L=101)   350     96        0 ( 582/ 747- 662/ 794)  18.4 E22-E16 S28-S24
     17 1667  507 (S16W60 L=180)   904    180        1 (1622/ 466-1696/ 539)  76.0 W55-W63 S18-S14
     18 1339  271 (S40W49 L=169)  1067    252        1 (1256/ 257-1437/ 288)  14.9 W42-W57 S42-S37
  

    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