Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  675 1518 (N33E28 L= 13)   690    133        0 ( 621/1469- 725/1561) -26.6 E33-E23 N29-N35
      2 1202 1329 (N22W11 L= 52)  1754    402        2 (1011/1226-1389/1435) -25.6 E02-W23 N16-N28
      3  307 1351 (N19E55 L=346)   902    177        0 ( 283/1289- 329/1430) -66.0 E57-E51 N15-N24
      4  717 1397 (N25E23 L= 18)   407     55        1 ( 692/1373- 749/1424) -26.6 E24-E20 N23-N26
      5  876 1299 (N19E11 L= 30)   769    232        1 ( 828/1212- 957/1406) -56.3 E14-E04 N14-N26
      6  729 1270 (N17E20 L= 21)   310     77        2 ( 694/1239- 757/1297) -53.1 E22-E17 N14-N18
      7  491 1140 (N08E35 L=  6)   935    207        0 ( 461/1047- 549/1232) -57.5 E38-E30 N02-N13
      8  403  975 (S03E41 L=360)   398    111        0 ( 366/ 955- 451/1006)  -8.1 E44-E37 S04-S01
      9 1011  499 (S31W01 L= 42)   319     74        0 ( 982/ 468-1035/ 532) -63.4 W00-W03 S33-S29
     10 1194  261 (S51W21 L= 62)   656    129        1 (1144/ 226-1230/ 290) -28.6 W15-W26 S54-S48
  

    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