Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1243 1337 (N16W04 L=185)   428    211        0 (1183/1257-1279/1415) -45.0 E01-W06 N11-N19
      2  219 1241 (S13E61 L=121)   455    195        0 ( 185/1210- 260/1289) -18.4 E67-E55 S14-S10
      3 1071 1195 (N04E02 L=180)   238     69        4 (1042/1161-1092/1221) -45.0 E04-W00 N01-N05
      4  352  938 (S28E42 L=140)   860    238        0 ( 276/ 879- 424/1007)  45.0 E47-E36 S33-S22
      5 1524  772 (S06W36 L=218)   312    106        2 (1476/ 747-1571/ 791)   8.1 W32-W39 S07-S04
      6  872  602 (S35W04 L=185)   267     53        2 ( 842/ 583- 892/ 627)   0.0 W01-W05 S36-S32
      7 1173  564 (S28W23 L=205)   651    229        0 (1124/ 489-1245/ 617) -77.5 W19-W28 S31-S22
      8 1645  387 (S21W66 L=248)  1857    285        4 (1585/ 352-1731/ 430)  28.1 W60-W75 S23-S15
  

    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