Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1376 1494 (N28W18 L=319)   321     67        0 (1345/1457-1401/1519) -53.1 W15-W19 N25-N29
      2  631 1381 (N12E28 L=272)  1441    290        4 ( 554/1297- 712/1503) -40.6 E36-E22 N07-N19
      3 1801 1156 (N13W52 L=352)  2261    417        5 (1712/1049-1912/1218) -15.9 W43-W65 N07-N16
      4 1437 1112 (N04W24 L=325)   323    109        2 (1389/1069-1483/1138)  -8.1 W20-W27 N02-N05
      5 1337 1023 (S02W19 L=319)   232     88        0 (1308/1000-1378/1063) -21.8 W16-W21 S03-N00
      6  166  923 (S17E62 L=238)  1850    276        3 ( 115/ 840- 211/1004) -39.8 E69-E57 S22-S12
      7  667  895 (S17E21 L=280)   283     70        1 ( 624/ 882- 711/ 911)  21.8 E23-E18 S17-S15
      8 1717  775 (S11W50 L=350)   524     86        1 (1687/ 738-1752/ 810) -23.2 W46-W53 S12-S09
      9  205  756 (S28E59 L=242)   407    117        1 ( 176/ 738- 233/ 776)  15.9 E62-E55 S28-S26
     10 1350  616 (S26W27 L=327)  1962    380        1 (1177/ 539-1547/ 748)  39.0 W16-W37 S32-S15
     11 1039  512 (S37W08 L=308)   448     94        2 ( 990/ 496-1092/ 538)  15.9 W04-W11 S38-S35
  

    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