Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1681 (N43E29 L=280)  3412    486        3 ( 546/1602- 970/1760)  -8.1 E59-E09 N38-N47
      2 1034 1551 (N36E04 L=306)   569     98        1 (1000/1514-1085/1580)  31.0 E06-W00 N32-N37
      3 1318 1531 (N37W18 L=327)  1411    269        4 (1185/1503-1435/1558)   2.9 W07-W28 N33-N38
      4 1564 1431 (N31W38 L=347)   926    163        4 (1506/1372-1627/1480) -33.7 W33-W42 N27-N34
      5  552 1380 (N20E35 L=274)   707    110        1 ( 522/1329- 586/1436) -74.1 E37-E32 N16-N23
      6  314 1380 (N18E56 L=253)  1173    123        3 ( 257/1342- 362/1424)  21.8 E62-E51 N15-N20
      7 1296 1309 (N21W16 L=325)   343    138        0 (1269/1263-1339/1374) -60.3 W13-W18 N18-N25
      8 1383 1318 (N22W22 L=331)   303     91        0 (1367/1271-1396/1369)  80.5 W20-W22 N19-N25
      9 1819 1247 (N20W60 L=  9)   398     37        1 (1801/1226-1836/1275) -56.3 W57-W62 N18-N21
     10  502 1120 (N04E34 L=275)  1491    250        2 ( 467/ 999- 549/1246)  85.9 E36-E30 S03-N11
     11  157 1209 (N06E71 L=239)   861     82        1 ( 130/1185- 196/1235)   9.5 E76-E64 N04-N07
     12  371 1000 (S04E44 L=266)   372     68        2 ( 337/ 981- 411/1018)   9.5 E46-E40 S05-S02
     13 1810  922 (N00W57 L=  7)   500     82        0 (1778/ 903-1849/ 951)   8.1 W54-W61 S01-N01
     14  207  803 (S18E62 L=247)   424     38        1 ( 193/ 778- 224/ 833)  56.3 E64-E59 S19-S16
     15  730  741 (S18E17 L=292)   410    172        1 ( 651/ 710- 790/ 784) -16.7 E22-E12 S19-S15
     16  531  671 (S24E32 L=278)   531    110        1 ( 480/ 639- 583/ 709)  36.9 E36-E27 S26-S20
     17 1492  627 (S20W35 L=344)  1958    336        3 (1354/ 554-1649/ 707)  28.8 W26-W46 S25-S13
     18 1200  498 (S30W17 L=326)   771    122        2 (1141/ 475-1260/ 523)  20.6 W12-W20 S32-S28
  

    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