Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  696 1470 (N30E25 L= 15)  1076    261        1 ( 622/1374- 754/1562) -40.2 E33-E19 N23-N35
      2 1217 1326 (N22W12 L= 52)  1636    380        2 (1041/1226-1390/1434) -27.6 W00-W23 N16-N28
      3  310 1348 (N19E55 L=346)   817    147        0 ( 287/1286- 331/1413) -63.4 E57-E51 N15-N23
      4  888 1289 (N19E10 L= 31)   703    224        0 ( 837/1212- 969/1388) -48.0 E13-E04 N14-N25
      5  735 1270 (N17E20 L= 21)   333     91        2 ( 697/1234- 763/1298) -31.0 E22-E17 N14-N18
      6  499 1137 (N07E35 L=  6)   901    215        0 ( 466/1047- 555/1234) -57.5 E37-E30 N02-N13
      7  409  972 (S03E41 L=360)   359    100        1 ( 373/ 953- 455/ 996)  -9.5 E43-E37 S04-S01
      8 1019  498 (S31W02 L= 42)   315     66        0 ( 992/ 468-1042/ 527) -53.1 W00-W03 S33-S29
      9 1199  260 (S51W21 L= 62)   649    126        0 (1147/ 225-1236/ 289) -26.6 W15-W27 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