Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  353 1436 (N01E55 L=357)   612     87        1 ( 317/1413- 398/1462)  -5.7 E60-E50 N00-N02
      2 1433 1402 (N26W14 L= 67)   277     68        0 (1399/1376-1463/1426) -11.3 W11-W16 N24-N26
      3 1222 1335 (N15W03 L= 55)   477    160        1 (1170/1272-1250/1398) -39.8 E01-W05 N12-N17
      4 1420 1145 (N10W19 L= 72)   210    138        0 (1398/1077-1437/1201) -50.2 W16-W21 N06-N12
      5 1793 1121 (N21W47 L=100)   411     64        0 (1777/1077-1804/1157) -53.1 W45-W48 N18-N23
      6  364 1093 (S18E44 L=  9)   970    205        0 ( 282/1058- 457/1127)  15.9 E52-E37 S20-S14
      7  273  853 (S35E51 L=  2)   668    118        0 ( 225/ 816- 317/ 885)  -5.2 E56-E45 S36-S33
      8  731  721 (S32E10 L= 43)   225     86        0 ( 711/ 680- 743/ 763) -59.0 E11-E08 S34-S29
      9  968  640 (S30W08 L= 60)   375    151        0 ( 902/ 620-1043/ 660)  16.7 W03-W13 S31-S28
     10  379  412 (S65E55 L=358)  1325    137        1 ( 312/ 389- 446/ 440)   2.5 E85-E38 S66-S63
     11  745  241 (S67W13 L= 65)   367     44        0 ( 727/ 219- 759/ 268) -35.5 W10-W17 S69-S64
     12 1062  178 (S57W50 L=102)   398     68        0 (1046/ 156-1088/ 205)  -9.5 W44-W56 S58-S54
  

    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