Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1542 1549 (N20W50 L=248)   546    120        1 (1501/1515-1594/1569) -33.7 W47-W54 N17-N21
      2  774 1382 (N33E05 L=193)   229     63        0 ( 747/1352- 794/1422)  31.0 E07-E02 N31-N35
      3  540 1306 (N35E24 L=174)   544    136        1 ( 480/1264- 611/1348)   5.2 E28-E17 N33-N36
      4  960 1182 (N17W01 L=198)   239     76        2 ( 925/1145-1000/1203) -78.7 E01-W01 N13-N18
      5  413 1117 (N27E37 L=161)   315     53        1 ( 382/1098- 439/1135)   0.0 E39-E34 N25-N27
      6  362  950 (N18E44 L=154)   302     60        0 ( 338/ 924- 390/ 980) -45.0 E45-E41 N16-N20
      7 1018  523 (S22E15 L=183)   267     65        0 ( 991/ 497-1053/ 552) -76.0 E15-E12 S24-S19
      8  492  440 (S14E52 L=146)  1668    319        4 ( 413/ 356- 539/ 530) -74.7 E57-E45 S20-S09
      9 1715  470 (S53W40 L=237)   616    114        1 (1679/ 456-1743/ 491) -24.4 W34-W45 S56-S49
  

    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