Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  565 1662 (N19E50 L=225)   389     41        0 ( 542/1641- 592/1688)  76.0 E51-E48 N16-N20
      2  351 1421 (S00E55 L=219)   424     66        0 ( 336/1387- 365/1451) -36.9 E57-E53 S02-N01
      3  813 1380 (N08E22 L=253)   593    174        0 ( 751/1323- 893/1448) -14.0 E27-E15 N06-N10
      4 1283 1229 (N12W09 L=283)  1233    321        1 (1245/1060-1306/1377) -60.9 W03-W13 N02-N20
      5  197 1193 (S16E62 L=213)   598    101        2 ( 173/1163- 224/1226) -63.4 E65-E59 S18-S14
      6 1493  952 (N03W28 L=303)   372     98        0 (1465/ 905-1513/1009) -71.6 W26-W29 N00-N06
      7  334  926 (S29E43 L=232)  1077    274        1 ( 237/ 892- 421/ 967)  14.0 E52-E36 S32-S25
      8  978  891 (S15W01 L=276)  1191    370        0 ( 791/ 851-1153/ 930)  15.9 E10-W11 S18-S11
  

    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