Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1866 (N41E51 L=344)  1356    154        2 ( 781/1829- 898/1894)   0.0 E58-E43 N37-N43
      2  439 1624 (N13E60 L=335)  1837    259        2 ( 390/1519- 489/1750) -54.8 E69-E55 N06-N23
      3 1231 1591 (N33E05 L= 30)   643    127        1 (1204/1537-1253/1655) -49.4 E09-E02 N29-N37
      4  358 1615 (N12E69 L=326)   658     67        1 ( 337/1586- 382/1653) -35.5 E72-E64 N09-N14
      5  475 1522 (N08E50 L=345)  2181    340        2 ( 364/1438- 605/1616)  51.8 E56-E44 N01-N16
      6  255 1152 (S18E56 L=339)   829    129        0 ( 193/1119- 305/1183)  22.6 E63-E51 S20-S15
      7  120 1014 (S28E75 L=320)   717     38        1 ( 100/ 996- 143/1031)  10.3 E80-E69 S29-S27
      8  588  849 (S28E22 L= 13)   902    208        5 ( 499/ 803- 685/ 906)  42.3 E27-E16 S32-S22
      9 1016  643 (S27W11 L= 46)   559    277        2 ( 888/ 615-1145/ 667)  35.5 W04-W18 S32-S22
     10 1391  603 (S18W35 L= 70)   990    189        1 (1310/ 558-1478/ 646)  59.0 W31-W38 S22-S12
     11  992  492 (S37W16 L= 51)  1156    272        1 ( 894/ 405-1095/ 556)  75.1 W12-W20 S45-S30
     12 1220  180 (S48W62 L= 97)   556     53        2 (1197/ 149-1244/ 208) -20.0 W57-W69 S49-S45
  

    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