Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1271 1680 (N35W38 L=141)   651    199        2 (1176/1648-1356/1709) -41.6 W32-W42 N30-N39
      2 1692 1636 (N19W70 L=174)  1417    118        3 (1666/1575-1718/1685)  11.9 W63-W83 N16-N21
      3 1318 1567 (N27W34 L=137)   693    181        0 (1255/1498-1387/1600)   9.5 W27-W39 N24-N28
      4 1402 1418 (N16W33 L=136)   249     82        1 (1369/1384-1419/1455) -90.0 W31-W34 N13-N18
      5 1131 1354 (N19W15 L=118)   193     85        0 (1099/1319-1157/1390) -78.7 W13-W15 N17-N22
      6  855 1286 (N23E03 L=100)   782    203        0 ( 803/1201- 937/1364) -71.6 E05-W00 N15-N27
      7 1745 1159 (S07W48 L=152)  2164    402        1 (1691/ 972-1784/1316)  63.4 W41-W53 S16-N02
      8  600 1217 (N25E21 L= 82)  2086    578        4 ( 367/1065- 796/1295) -35.1 E39-E10 N11-N32
      9  491 1081 (N19E31 L= 72)   227    109        0 ( 437/1066- 542/1097) -29.7 E34-E27 N17-N21
     10 1552  847 (S19W27 L=130)  1553    316        4 (1425/ 746-1679/ 920)   8.1 W16-W37 S21-S16
     11 1842  861 (S27W53 L=157)  2354    227        5 (1733/ 805-1932/ 920)  -4.4 W41-W68 S30-S24
     12  368  600 (S05E54 L= 50)  5599    464       13 ( 254/ 446- 463/ 830)  55.7 E64-E42 S15-N07
     13  757  718 (S07E23 L= 81)   325    107        0 ( 724/ 674- 776/ 766)  45.0 E25-E20 S08-S04
     14 1613  715 (S29W30 L=133)   325     73        0 (1576/ 688-1638/ 750)  18.4 W26-W32 S30-S27
     15  339  644 (S02E54 L= 49)   375     72        1 ( 315/ 607- 364/ 678) -76.0 E55-E52 S04-N00
     16  610  328 (S27E50 L= 53)   416     59        3 ( 578/ 308- 647/ 346)  -8.1 E53-E46 S28-S25
     17 1414  231 (S63W05 L=109)   453     46        0 (1386/ 215-1442/ 250) -45.0 W02-W09 S66-S59
  

    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