Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  551 1290 (N22E28 L=309)  1615    325        5 ( 446/1189- 651/1368) -29.4 E37-E21 N14-N28
      2  127 1205 (N22E72 L=265)   880    204        1 (  91/1171- 164/1230)   7.1 E86-E64 N20-N25
      3  283 1186 (N19E50 L=287)   299    139        1 ( 239/1170- 327/1201) -11.3 E55-E46 N17-N21
      4 1224 1058 (S01W12 L=349)   302    147        0 (1156/1016-1297/1083)   6.3 W08-W16 S04-N01
      5 1366  994 (S06W20 L=357)   310    134        1 (1315/ 962-1429/1022) -33.7 W17-W24 S09-S04
      6 1309  948 (S08W16 L=353)   312    135        0 (1248/ 934-1378/ 969) -20.6 W13-W20 S10-S06
      7 1182  866 (S12W07 L=344)   349    121        0 (1147/ 801-1215/ 925)  50.2 W06-W09 S16-S08
      8 1024  798 (S14E03 L=334)   558    138        0 (1009/ 728-1039/ 873)  69.4 E03-E03 S18-S09
      9 1587  660 (S30W36 L= 13)  1609    324        4 (1491/ 547-1639/ 767)  65.2 W27-W44 S39-S22
     10  163  704 (S08E74 L=263)  2026    208        0 ( 137/ 642- 187/ 753)  45.0 E75-E74 S12-S05
     11  529  568 (S21E41 L=296)  1165    237        1 ( 482/ 482- 566/ 671)  37.6 E41-E42 S27-S14
     12  398  541 (S21E54 L=283)   309     63        0 ( 368/ 517- 431/ 568) -45.0 E55-E52 S23-S19
  

    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