Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1122 1838 (N53E26 L=154)  1607    187        1 (1037/1799-1191/1867)  17.1 E32-E17 N48-N56
      2  824 1793 (N38E42 L=137)   377     26        0 ( 803/1775- 847/1808)  90.0 E44-E40 N36-N39
      3 1317 1703 (N48E01 L=178)   378     57        0 (1288/1682-1346/1722)  53.1 E03-W00 N45-N49
      4  458 1617 (N17E58 L=122)   829    117        2 ( 432/1587- 506/1665) -15.3 E63-E52 N14-N19
      5  500 1503 (N11E47 L=133)   434     61        1 ( 474/1476- 527/1534)  78.7 E48-E45 N08-N13
      6 1747 1409 (N41W45 L=224)   516     84        1 (1698/1397-1786/1423)  18.4 W38-W50 N38-N43
      7 1144 1398 (N20E03 L=177)   429    153        1 (1062/1366-1206/1415)  14.0 E07-W01 N18-N21
      8 1222 1202 (N11W07 L=186)   200     57        1 (1188/1188-1251/1218)  36.9 W04-W08 N09-N12
      9 1748 1065 (N18W46 L=226)   546     64        1 (1722/1025-1765/1101) -78.7 W43-W47 N15-N20
     10  393  997 (S21E40 L=140)   469    169        0 ( 369/ 932- 408/1081) -71.6 E41-E38 S24-S15
     11  912  930 (S12E04 L=176)   552    194        2 ( 874/ 837- 937/1015) -85.2 E05-E02 S18-S06
     12  516  916 (S23E29 L=151)   316     61        1 ( 482/ 890- 540/ 940)  14.0 E31-E27 S24-S21
     13 1181  815 (S12W14 L=194)   268     51        0 (1154/ 794-1210/ 835)   0.0 W12-W16 S12-S10
     14  436  734 (S36E33 L=146)   485    108        1 ( 394/ 698- 471/ 772)  45.0 E37-E30 S39-S32
     15  986  661 (S26W08 L=187)   249     86        0 ( 943/ 641-1029/ 682)   0.0 W04-W10 S26-S25
     16 1663  492 (S16W62 L=242)   495     57        0 (1634/ 466-1688/ 519)  76.0 W60-W64 S17-S13
     17 1490  423 (S25W52 L=231)   666    146        0 (1432/ 373-1546/ 477)  83.7 W50-W53 S29-S20
  

    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