Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1647 1553 (N23W55 L=221)   834    156        0 (1588/1494-1706/1600) -58.0 W52-W58 N18-N26
      2  626 1510 (N36E21 L=145)   793    232        0 ( 544/1450- 736/1543)   3.8 E27-E12 N33-N38
      3 1264 1101 (N01W15 L=181)   192     92        0 (1239/1067-1298/1144)  36.9 W13-W17 N00-N03
      4 1561  949 (S12W32 L=198)   269    129        0 (1536/ 886-1580/1003)  60.3 W30-W34 S15-S08
      5 1196  724 (S20W06 L=172)   554    216        0 (1104/ 661-1254/ 799)  31.0 W00-W10 S22-S16
      6 1487  746 (S23W26 L=191)   221     53        0 (1475/ 717-1500/ 782) -76.0 W25-W26 S24-S20
      7  993  661 (S21E07 L=158)   242     86        0 ( 953/ 644-1040/ 681) -31.0 E09-E04 S22-S19
      8  922  636 (S21E12 L=153)   251     64        0 ( 909/ 595- 938/ 668)  53.1 E13-E10 S23-S19
      9  347  609 (S15E54 L=111)   583    106        0 ( 331/ 561- 371/ 665)  80.5 E56-E52 S17-S11
     10  292  607 (S14E60 L=105)  1161     95        2 ( 257/ 558- 318/ 658)  60.3 E64-E56 S17-S10
     11  686  599 (S20E28 L=137)   302     76        0 ( 662/ 561- 714/ 635)  38.7 E30-E25 S22-S18
     12  418  375 (S31E64 L=102)  1706    248        1 ( 347/ 319- 509/ 456) -51.3 E70-E55 S35-S24
  

    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