Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  958 1827 (N64E08 L=  5)   605     77        1 ( 927/1796- 996/1846)  18.4 E11-E02 N60-N66
      2  969 1600 (N42E04 L=  9)   717    174        0 ( 910/1535-1022/1672) -50.7 E09-W00 N37-N48
      3 1509 1504 (N34W39 L= 52)   283    109        0 (1479/1457-1540/1540) -71.6 W37-W40 N30-N36
      4 1183 1329 (N23W11 L= 24)   303    150        0 (1128/1282-1235/1373) -40.6 W07-W14 N19-N25
      5  205 1317 (N20E69 L=304)   615     83        0 ( 179/1292- 243/1348)  -7.6 E77-E62 N19-N21
      6 1351 1317 (N22W22 L= 35)   238     84        0 (1307/1305-1383/1339)  11.3 W19-W24 N20-N22
      7  530 1268 (N19E34 L=339)   247     66        1 ( 498/1250- 573/1282)  -9.5 E36-E30 N17-N19
      8 1418 1146 (N11W26 L= 39)   216    102        0 (1391/1102-1441/1192) -71.6 W24-W26 N07-N13
      9 1700  720 (S17W49 L= 62)  3735    597        7 (1597/ 553-1802/ 892)  74.7 W42-W57 S28-S06
     10 1611  759 (S14W40 L= 53)   831    101        1 (1567/ 726-1661/ 805)  15.9 W36-W44 S16-S11
     11 1231  696 (S17W13 L= 26)   703    201        1 (1150/ 642-1316/ 743)  31.0 W08-W18 S20-S14
     12  646  493 (S31E29 L=344)  2385    300        5 ( 567/ 404- 726/ 590) -56.3 E34-E22 S38-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