Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  852 1465 (N31E13 L= 38)   377     83        1 ( 820/1437- 897/1483)  18.4 E15-E09 N29-N32
      2 1541 1319 (N22W36 L= 87)   302    108        0 (1497/1286-1573/1340) -38.7 W32-W38 N19-N23
      3 1140 1250 (N17W07 L= 58)   387     96        0 (1092/1231-1180/1269)  11.3 W04-W09 N16-N18
      4  420 1207 (N13E42 L=  9)   482    144        0 ( 392/1138- 457/1265) -48.8 E45-E38 N08-N16
      5 1771  853 (S08W54 L=105)   412     79        0 (1746/ 821-1801/ 880)  45.0 W51-W57 S09-S05
      6  243  825 (S12E58 L=353)   899     97        2 ( 201/ 789- 276/ 850) -20.6 E63-E54 S13-S10
      7  612  559 (S28E29 L= 22)   256     61        0 ( 602/ 526- 623/ 594)  78.7 E30-E28 S30-S25
      8 1435  422 (S37W35 L= 86)   981    201        1 (1406/ 337-1460/ 516) -68.2 W31-W40 S44-S29
  

    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