Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1919 (N69W17 L=331)  1435    139        3 ( 889/1871- 968/1957)  24.6 W08-W34 N62-N75
      2 1322 1599 (N31W30 L=345)   376    106        0 (1276/1560-1360/1629) -45.0 W27-W31 N28-N33
      3 1349 1435 (N20W27 L=342)   757    266        0 (1235/1346-1426/1504) -55.3 W20-W29 N13-N26
      4  849 1482 (N30E05 L=310)   385     98        0 ( 797/1464- 900/1498) -23.2 E08-E01 N28-N31
      5  461 1275 (N22E34 L=281)   488     74        1 ( 440/1238- 483/1319)  53.1 E35-E32 N19-N25
      6  232 1207 (N21E55 L=260)   386     56        2 ( 216/1177- 253/1232)  56.3 E57-E53 N19-N22
      7 1363 1115 (N01W21 L=336)   269     90        0 (1325/1078-1397/1145)  21.8 W18-W23 N00-N02
      8  588  576 (S21E36 L=279)  1882    312        4 ( 541/ 448- 665/ 726)  78.1 E38-E32 S30-S11
      9 1648  572 (S37W42 L=357)  2863    434        4 (1481/ 472-1799/ 698)  19.0 W26-W56 S41-S30
     10 1027  479 (S34E08 L=307)   353    151        1 ( 943/ 465-1089/ 495)  -5.2 E14-E03 S35-S32
     11  713  394 (S35E34 L=281)   314     80        0 ( 699/ 356- 736/ 424)  33.7 E37-E30 S37-S33
     12  927  282 (S47E23 L=292)   307     83        0 ( 892/ 269- 976/ 296)   0.0 E26-E18 S48-S45
  

    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