Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1464 1691 (N29W55 L= 81)   558     86        1 (1431/1671-1512/1714) -51.3 W50-W58 N26-N31
      2  548 1563 (N51E20 L=  7)   399     29        1 ( 520/1542- 577/1581) -18.4 E22-E16 N49-N53
      3 1139 1466 (N26W20 L= 46)   209     66        0 (1106/1454-1175/1478) -21.8 W17-W22 N24-N27
      4  212 1311 (N42E60 L=326)   712    104        0 ( 162/1284- 257/1327) -11.9 E71-E52 N39-N43
      5 1057 1234 (N16W07 L= 34)   439    174        0 ( 973/1206-1136/1257) -36.9 W03-W11 N12-N18
      6  545 1165 (N25E26 L=  1)   349     86        0 ( 506/1136- 590/1183) -33.7 E28-E22 N22-N26
      7  185 1110 (N30E61 L=325)  1224    175        0 ( 124/1034- 245/1166) -27.8 E72-E53 N24-N34
      8  345 1002 (N20E44 L=342)   434     82        0 ( 298/ 983- 392/1021)  -6.3 E48-E39 N19-N20
      9 1437  515 (S37W12 L= 38)   934    162        1 (1381/ 450-1486/ 587)  33.7 W06-W16 S41-S32
     10  471  489 (S12E52 L=335)   474     66        2 ( 441/ 459- 498/ 517) -53.1 E53-E49 S14-S10
     11  878  345 (S31E31 L=355)  2120    385        2 ( 738/ 220- 972/ 482)  19.5 E50-E19 S36-S23
     12 1203  375 (S39E09 L= 17)  1261    282        1 (1078/ 322-1322/ 435)  -2.9 E19-W01 S41-S36
  

    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