Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1668 1653 (N30W74 L= 54)   698     70        0 (1626/1636-1703/1670) -16.7 W68-W79 N28-N31
      2 1399 1549 (N31W38 L= 18)   456     64        0 (1369/1523-1436/1568)   0.0 W34-W41 N29-N32
      3  674 1335 (N31E19 L=322)  1743    464        5 ( 486/1250- 870/1442) -31.0 E35-E05 N22-N40
      4 1130 1133 (N11W08 L=349)   273     59        1 (1108/1100-1152/1168)  53.1 W06-W09 N09-N13
      5  456 1042 (N15E37 L=304)   332    134        0 ( 420/ 995- 513/1072)  29.7 E39-E32 N12-N16
      6  885 1021 (N09E08 L=332)   245     80        0 ( 867/ 975- 904/1063)  63.4 E09-E07 N06-N10
      7  268  980 (N13E54 L=287)  1011    197        2 ( 206/ 908- 327/1048) -36.9 E60-E48 N08-N17
      8 1190  801 (S09W06 L=347)   637    146        0 (1170/ 736-1211/ 873)  77.5 W04-W08 S13-S04
      9 1410  734 (S17W19 L=360)   485    174        1 (1372/ 695-1470/ 771) -26.6 W16-W23 S19-S14
     10  666  599 (S14E30 L=311)   782    302        1 ( 615/ 544- 737/ 661) -76.0 E32-E24 S17-S09
     11  582  523 (S18E38 L=302)   292     64        2 ( 564/ 491- 601/ 547)  31.0 E40-E35 S20-S16
     12  886  394 (S31E21 L=319)  1092    229        1 ( 805/ 334- 973/ 483)  26.6 E29-E13 S35-S25
  

    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