Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  597 1739 (N43E42 L=184)   365     30        0 ( 575/1717- 618/1762)  56.3 E43-E40 N40-N44
      2 1460 1643 (N39W32 L=258)  2291    434        2 (1335/1531-1620/1716) -30.1 W21-W43 N31-N44
      3 1229 1573 (N32W12 L=238)   275     77        1 (1194/1539-1262/1604) -38.7 W09-W14 N29-N33
      4  346 1491 (N24E53 L=173)   975    165        4 ( 325/1414- 371/1566) -57.5 E59-E50 N18-N29
      5  937 1531 (N27E08 L=218)   220     64        0 ( 903/1513- 972/1548) -21.8 E10-E05 N26-N28
      6 1696 1422 (N24W48 L=274)   469     94        1 (1674/1383-1717/1473) -78.7 W44-W49 N22-N27
      7 1550 1179 (N08W33 L=259)   330     96        1 (1534/1133-1566/1229) -78.7 W31-W33 N05-N10
      8  425 1091 (S02E39 L=187)   235     69        2 ( 408/1054- 440/1126) -71.6 E40-E37 S03-N00
      9  269  901 (S13E53 L=173)   753    189        3 ( 250/ 826- 292/ 985)  78.7 E55-E50 S17-S07
     10  510  911 (S12E33 L=193)   365    125        0 ( 479/ 867- 547/ 966)  50.2 E35-E29 S14-S08
     11 1088  806 (S17W05 L=231)   221     33        1 (1069/ 787-1111/ 825) -33.7 W03-W06 S18-S15
     12  934  788 (S19E05 L=221)   248    123        2 ( 896/ 763- 986/ 824)  38.7 E07-E01 S20-S16
     13 1731  755 (S16W51 L=278)   562     97        3 (1705/ 718-1761/ 796)  68.2 W49-W54 S18-S13
     14  472  626 (S30E39 L=187)   955    157        3 ( 443/ 554- 511/ 708) -90.0 E41-E35 S34-S24
     15  215  609 (S30E71 L=155)  1171     80        2 ( 188/ 570- 247/ 648)  12.5 E83-E64 S32-S28
     16 1384  590 (S29W27 L=253)   546    132        1 (1336/ 559-1412/ 643)  59.0 W24-W29 S31-S25
     17 1143  590 (S30W10 L=236)   575    221        1 (1038/ 572-1244/ 629)  12.1 W02-W16 S31-S27
     18  879  496 (S38E08 L=218)   266     65        0 ( 844/ 473- 906/ 524)  31.0 E11-E06 S39-S35
     19  616  302 (S54E40 L=186)   837    147        3 ( 551/ 256- 691/ 348)  15.1 E54-E28 S58-S50
  

    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