Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  318 1401 (S02E58 L=191)   465     74        0 ( 292/1366- 348/1435)  63.4 E59-E55 S04-N00
      2  643 1338 (N01E31 L=218)   436    140        0 ( 612/1270- 672/1391) -29.7 E33-E26 S02-N03
      3 1620 1080 (N14W33 L=282)  2255    283        4 (1557/ 885-1661/1339) -53.5 W23-W41 N04-N27
      4 1214 1188 (N08W06 L=255)   682    242        0 (1122/1109-1317/1283) -15.9 E01-W13 N05-N10
      5  143 1238 (S14E71 L=177)   761    120        2 ( 120/1205- 168/1273) -76.0 E74-E67 S16-S12
      6  365 1093 (S19E43 L=205)  1337    515        4 ( 202/1046- 561/1144)  15.4 E58-E29 S22-S13
      7  186  818 (S39E57 L=192)  1738    403        3 ( 125/ 755- 248/ 856)  28.3 E63-E50 S43-S36
      8  325  813 (S37E42 L=206)   328     74        0 ( 304/ 782- 340/ 855)  68.2 E44-E41 S38-S33
      9 1020  782 (S20W07 L=255)   241     77        1 ( 983/ 770-1066/ 799)  21.8 W04-W09 S20-S18
     10  529  606 (S45E22 L=227)   572    138        0 ( 473/ 588- 580/ 635)  33.7 E27-E17 S47-S41
  

    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