Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1423 1703 (N36W53 L=295)  1395    257        1 (1309/1667-1521/1727) -31.0 W43-W60 N30-N39
      2  409 1479 (N49E41 L=200)   502     87        0 ( 383/1450- 444/1517) -29.1 E45-E36 N46-N51
      3 1685 1415 (N12W55 L=296)   423     88        0 (1642/1378-1717/1440) -59.0 W51-W56 N09-N14
      4 1398 1286 (N13W29 L=270)  1017    245        1 (1361/1175-1444/1390)  81.9 W26-W32 N05-N19
      5  772 1076 (N15E14 L=228)  1183    515        0 ( 645/ 905- 884/1222) -73.1 E20-E09 N03-N26
      6 1022  894 (S00E03 L=238)  1142    381        0 ( 930/ 749-1156/ 951)  17.4 E11-W05 S06-N03
      7  137  907 (N15E72 L=169)   957     62        2 ( 112/ 880- 159/ 949) -76.0 E77-E68 N13-N17
      8 1829  822 (S26W53 L=295)   874     96        1 (1809/ 776-1851/ 865)  68.2 W50-W55 S29-S23
      9  376  552 (S09E56 L=186)  2687    495        6 ( 266/ 395- 452/ 696) -87.3 E63-E50 S20-N01
     10  298  526 (S10E67 L=175)  4746    367        9 ( 208/ 396- 390/ 663) -83.7 E73-E62 S19-N00
     11  524  582 (S09E42 L=199)   731    132        2 ( 493/ 533- 570/ 640) -66.8 E44-E39 S12-S04
     12 1679  528 (S42W38 L=279)  1712    269        4 (1566/ 449-1758/ 605)  14.7 W26-W47 S45-S37
     13  893  499 (S21E21 L=221)   611    126        1 ( 859/ 459- 922/ 567)  39.8 E23-E16 S24-S17
     14  766  487 (S19E29 L=212)   362     86        3 ( 729/ 453- 802/ 519)  15.9 E32-E25 S21-S17
     15  548  393 (S22E50 L=192)  1041    243        2 ( 485/ 306- 608/ 489)  22.9 E65-E39 S27-S16
     16 1264  203 (S57E08 L=233)   494     69        1 (1242/ 184-1290/ 238)  45.0 E10-E05 S60-S53
  

    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